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PREFACE 

The Agriculture and Resources Inventory Surveys Through Aerospace 
Remote Sensing progranii AgRlSTARS* Is a six-year program of researchi 
development, evaluation, and application of aerospace remote sensing for 
agricultural resources, which began In Fiscal Year l^oO. This program 
Is a cooperative effort of the National Aeronautics and Space Admini- 
stration. the U.S. Agency for International Development, and the U.S. 
Departments of Agriculture. Commerce, and the Interior. AgRlSTARS con- 
sists of eight Individual projects. 

The work reported herein was sponsored by the Supporting Research 
(SR) Project under the auspices of the National Aeronautics and Space 
Administration, NASA. Mr. Roberts. MacDonald, NASA Johnson Space Center, 
was the NASA Manager of the SR Project and Dr. Glen Houston was the 
Technical Coordinator for the reported effort. 

The Environmental Research Institute of Michigan and the Space 
Sciences Laboratory of the University of California at Berkeley comprise 
a consortliim having responsibility for development of corn/soybeans area 
estimation procedures applicable to South America within both the Sup- 
porting Research and Foreign Commodity Production Forecasting Projects 
of AgRlSTARS. 

This reported research, the extension of reflectance modeling to 
Include row effects, was performed within the Environmental Research 
Institute of Michigan's Infrared and Optics Division, headed by Richard 
R. Legault, Vice-President of ERIM, under the technical direction of 
Robert Horvath, Program Manager, and Dr. William A. Mallla, Task Leader. 
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1.0 

INTRODUCTION 

The scattering and reflection of Incident daylight by agricultural 
crops alters the direction and spectral composition of daylight radiation 
In a complex manner. A part of this altered radiation Is detected by 
Landsat. From the detected part of the altered radiation. It Is desired 
to infer the existence of important agronomic features of the crop which 
are the cause of the altered radiation. 

A canopy reflectance model provides the logical connection between 
the botanical features of the canopy, the geometry of the radiometric 
Interaction and the resulting alteration In the reflected radiation. 

Such a model allows one to understand the reasons for the alterations 
and to calculate the magnitude and trends of these alterations caused 
by the botanical features and the geometry of the Interaction. The 
validity of Inferences as to the existence of Important agronomic features 
from the detected altered radiation may be tested on theoretical grounds. 

Many crops are planted in rows by machinery. Upon emergence of the 
plants, the bare soil between rows is still the dominant feature which 
reflects incident daylight. As growth continues, the vegetation grows 
both higher and spreads out over the inter-row regions covering the bare 
soil. At some time during the grow^-ng season, the bare soil is covered 
enough that the bare soil between rows is no longer a dominant feature. 

The vegetation canopy becomes essentially laterally uniform in its 
radiation scattering properties. The alteration of incident daylight 
can be understood and calculated by a previously developed uniform canopy 
reflectance model [1] at this stage of growth. 

However, for a considerable time during the early part of the growing 
season, the strips of bare soil between rows and the increasing density 
of vegetation along the rows become equally Important in their contributions 
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to canopy reflectance. One may intuitively understand that the direction 
of sunlight relative to the row direction will change the relative influence 
of vegetation and bare soil. When the sun is directed along the row 
direction, the bare soil is fully illuminated but, when the sun is directed 
across rows, the soil is largely in the shadow of the standing vegetation 
along the rows. Thus, Landsat can receive different alterations due only 
to the way the rows trend relative to sunlight. An inference that such 
altered radiation is due to a change in some ioq>ortant agronomic feature 
could be in error. 

Field measurements of wheat [2] and soybeans [3] show that row 
direction relative to sunlight does change the character of the scattered 
and reflected radiation. The research presented in this report is the 
extension of the uniform canopy model so that it will also apply to the 
early season nonuniform row crop canopy so that this "row effect" may be 
understood and calculated. 

Verhoef and Bunnlk [4] developed a "row effect" model by assiming 
that the vegetation along rows formed a rectangular prism of plant material 
with bare soil between. The limitation of the model to a particular plant 
profile does not permit one to calculate Intermediate amounts of soil cover 
between rows as the crop continues to grow. It is particularly desirable 
that any lateral distribution of vegetation should be allowed so that the 
model will apply to any crop at any time in the growing season. 

The following text reviews the concepts, nomenclature, and symbols 
of the uniform canopy model in order to form the logical basis for its 
modification to incorporate the "row effect." The concept of density 
modulation is introduced to account for the row structure of a canopy 
and the manner of calculation using such a concept is described. 

A comparison of the row effects predicted by the Verhoef-Bunnik 
rectangular prism model and the extended uniform canopy model is made 
for the extreme case of the rectangular prism row to show that the results 
are comparable in this case. 
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Finally, thv extended model Is epplled to wheet where more realistic 
vegetation distributions are assumed. The results are similar to those 
of field meacureokents. The red band> Landsat band 5» Ir most sensitive 
to row direction becau&e of the usual large contrast between vegetation 
and soil. Reflectance In this band may easily vary by a factor of two 
with changing r w direction. The IR bands, Landsat 6 and 7, are least 
.'fccccd by row direction because of low contrast between soil and vegeta- 
tion because of the large amount of diffuse flux scattered to soil by 
Che veaeCfttion. 
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2,0 

REVIEW OF THE UNIFORM CANOPY MODEL 

The uniform canopy reflectance model consists of a number of infinitely 
extended horizontal layers or strata as illustrated in Figure 1. Within each 
layer » the plant components of the canopy are considered to be randomly 
distributed and homogeneously mixed. The plant components are the identi** 
fiable parts of the plant* such as* stems* leaves* branches* flowers and 
pods. Most vegetation has a morphology which is largely genetically con- 
trolled. As a consequence* agricultural vegetation tends to be stratified 
in such a way th^t particular plant components occur at different levels 
above the soil. For example, the heads of wheat are located at the top 
of the canopy and not at the bottom. The influence of the heads of wheat 
on the reflectance of a wheat field will be greater because the top layer 
receives full Illumination and may be seen without obstruction. Lower 
layers lie in partial shadow and must be seen through the obstruction of 
the upper layers. 

The location of layers are chosen in siich a way as to quantize the 
vertical distribution of components. Each layer should contain uniformly 
some fraction of component types that is characteristic for that level 
above the soil. 

Collimated radiation from the sun enters the top of the canopy. 

This ccllimated flow of radiation is called specular flux in the following 
text. That specular flux which is Intercepted by a plant component is 
diffusely scattered and partially absorbed. The remaining specular flux* 
steadily diminished by such scattering* proceeds on to the soil making 
"sun flecks" upon the soil surface. 

The diffuse flux created by scattering may be produced by reflection 
from a component or by transmission through a component* Some of the 
diffuse flux is scattered towards the top of the canopy; the remainder 
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FIGURE 1. SCHEMATIC OF UNIFORM CANOPY MODEL. 
The actual plant components are illustrated on 
the left and the model equivalent components 
are shown on the right. 
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i* scattered towards the soil. As the diffuse flux swves through the 
canopy, some of the diffuse flux will be intercepted and scattered again 
with some of the rescattered flux going up and some going down and so forth. 

The lateral average fl\ix density on a horisontal plane of specular 
flux, and upward and downward welling diffuae flux, varies with depth in 
the canopy. Allen, Gayle, and Richardson [5] showed by experiment that 
the flux densitiea could be derived using Duntley'a differential equations 
for scattering in diffuse optical media. The scattering propertler of 
any particular medium are specified by the values assigned to five indepen- 
dent parameters in these equations. These differential equations are 
shown in relations (1), (2), and (3), 


dE(4d)/dz « -aE(4d) 

«*• bE(-d) + cE(s) 

(1) 

dE(-d)/dz - aE(-d) 

- bE(4d) - c'E(») 

(2) 

dE(s)/dz - k(Es) 


(3) 


where E(-hl) * upward welling diffuse flux density, 

E(-d) » downward welling diffuse flux density, 

E(s) ■ specular flux density, 

a • extinction coefficient for diffuse flux, 
b * backscatterlng coefficient for diffuse flux, 
c • backscatterlng coefficient for specular flux, 
c' ■ forward scattering coefficient for specular flux, 
k * extinction coefficient for specular flux. 

The five parameters, a, b, c, c', and k for each layer plus the boundary 
conditions of soil reflection at the bottom and sunlight at the top are 
all that Is needed to specify how much flux goes which way. What remains 
unknown Is the relationship between these parameters and the plant compo- 
nents that are present within the -anopy. 
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The uniform canopy modal provides a systematic end logical method 
of cslculeting epproxlMte values for these parasMters given the nui^r, 
orientation* and spectral properties of the plant components In a canopy. 
This method conccpf.ually replaces a particular plant component with 
three plane orthogotoal projections of that component. Each plane pro- 
jection (hereafter called a 0K>del equivalent cooq>onent) Is assigned the 
same hemispherical apectral reflectance and transmittance as that of 
the actual plant coo^onent. The concept of projections Is illustrated 
In Figure 2. The projection on the horizontal plane becomes the hori- 
zontal owdel equivalent component and the projections on vertical planes 
become the vertical model equivalent components. If the plant components 
were replaced by these model equivalent components* then collimated flux 
Incident on model equivalent cosq>onents in any of these three orthogonal 
directions would scatter flux In the same amount as would be scattered 
by the plant component Itself. For other Incident directions* the 
amount would only be approximately the same. 

The five unknown parameters can now be calculated using model equi- 
valent components. For one type of plant component* 


a 

b 


c 


c 

k 


h"h‘l 

- t) + “ Ip + t1/2)* 

('0 

hV 

+ o^n^(p/2 + t/2)* 

. \ 
/ 

'hV 

+ (2/it) u^«^(p/2 + t/2) tan 6^* 

(6) 

hV 

+ (2/it) o^n^(p/2 + t/2) tan 0^* 

(7) 


(2/tt) 0 n tan 6 * 
v V a 

(8) 


where o. - mean horizontal model equivalent component area* 
n 

0 ^ * mean vertical model equivalent component area* 
n^ - nusd>er of horizontal coisponents per unit volume* 
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FIGURE 2. THE USE OF PROJECTIONS AS MODEL EQUIVALENT COMPONENTS. 
Three orthogonal projections of a leaf show the relationship 
between a plant component (leaf) and the model equivalent. 
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- number of vertical components per unit volume • 

T ■ hemispherical transmittance of the plant component, 
p ■ hemispherical reflectance of the plant component, 

6 ■ polar angle of incident specular flux. 

s 

The mean values are taken over a field area the siee of the instantaneous 
field of view (IFOV) . The parameters for each type of plant component 
are determined in a similar manner. The parameter value for the ensemble 
of components within any given layer is simply the sum of the parameters 
for each type of component in that layer. Thus, for example, 

a(layer) ■ a(type 1) + a(type 2) + a(type 3) + .... 

and similarly for the other parameters. 

The relations (4) through (8) show how each plant component Influences 
the radiant flux within the canopy. The influence Increases as the size 
and concentration Increases. If the plant component changes orientation 
with the vertical direction, the model equivalent components, which are 
projections, change size. The influence of plant morphology is accounted 
for in this simple manner. 

Equipped with the values for the five parameters for each layer, one 
may solve relations (1), (2), and (3) for each layer and, hence, for the 
flux within the canopy. This flux is the llluminant for objects within 
the canopy which one can see from some direction of view. The final 
computation now is simply to determine the radiance, L, (radiometric 
brightness) of each component in the canopy and what fraction of these 
components can be seen without obstruction. The model equivalent compo- 
nents are again used to calculate the expected radiance of the components. 


10 




At any particular level within the canopy, the radiance per unit depth 

* 

is 

dL/dz ■ (uE(4d) + vE(-d) + w'E(s))/it, (9) 

where u ■ ^/2)(2/t) tan 6^, 

^ " ‘’hV’ ^/2)(2 /it) tan 0^, 

w' ■ %”v % C(sin 4* + cos <P)P 

+ (sin il* * Ip cos >i')T3/2Tr, 

<1' ■ azimuthal angle between sun and view directions, 

■ polar view direction. 

V 

The fraction , *) of components that can be seen at depth, z, 

from view angle, is given by 


dPv('^>v^'^‘** * *^^v^V ’ 

where K ■ Oj^n^^ + (2/i>) o^n^ tan 0^. 

Adding all contributions to radiance through the canopy and adding 
the soil radiance at the bottom of the canopy (z “ -hX all in proportion 
to the fraction within clc-"* view, one obtains the radiance that can be 
observed from the position of view, 

o 

L ■ X p (0 , z)[]uE(-Hl) + vE(-d) + w*E(s)] dz/n 

( 11 ) 

+ P^(%. -h)[p (soll)(E(-d) + E(s))]/w. 


The relation for w" here differs from that in the original publication 
which was in error. The difference is due to a missing factor which was 
pointed out by Verhoef and Bunnik and confirmed by W. Nalila. 
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The canopy reflectance for sunlight alone Is, 

p(sun, canopy) • irL/E(s). (12) 

The reflectance of the canopy for skylight alone is calculated in 
the same manner except that E(s) - 0 and the value of the downward welling 
diffuse flux, E(-d), at the top (z « 0) of the canopy, instead of being 
set to zero, is set equal to the value of downward welling diffuse sky- 
light, E(sky). 

The canopy reflectance for skylight alone is, 

p(sky, canopy) ■ xL/E(sky). (13) 

When both skylight and sunlight are incident upon the canopy together 
to produce global illumination, E(global), then the reflectance for global 
daylight is the weighted average of relations (12) and (13), 

p (global, canopy) ■ [E(sky)p (sky, canopy) +E(s)p(sun, canopy) ] /E (global) (14) 
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3.0 

EXTENSION TO INCLUDE ROW EFFECTS 



The fundamental concepts* nomenclature* and procedures of the uniform 
canopy model will be used with certain modifications to incorporate the 
effects of a row structure in agricultural crops. These modifications are 
introduced in such a way as to reduce to the uniform canopy model as row 
structure disappears from the crop due to overgrowth of the area between 
rows by the natural growth of the crop during the growing season. 

3.1 THE CONCEPT OF DENSITY MODULATION 

In the uniform canopy* the density of components are the mean values 
for a patch of field the size of the instantaneous field of view (IFOV). 
Locally* the densities can be expected to vary due to the randomness of 
the distribution. Random distributions are expected to be clumpy but 
without any order as to where the clumps occur. One could consider any 
narrow strip of field and determine the mean density of components within 
that strip. The mean density would be the same as the IFOV mean, given 
sufficient strip length for any direction the strip might take over a 
uniform canopy. 

However, in the case of a canopy with row structure, the strip mean 
will converge to a different mean density for strips parallel to the row 
direction depending upon the lateral displacement of the strip from the 
row center. The variation of strip means would be periodic for displacements 
of the strip in the across-row direction with large values on the row 
centers and small values between row centers. This variation in strip 
means relative to the IFOV mean is hereafter called density modulation. 
Density modulation is the evidence for the existence of row structure and 
is the measure of the amount of row structure. 

To describe crop row structure by geometric shapes, such as* rectan- 
gular prisms, semicircular prisms, or other specialized prism shapes, is 
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Co utilize special cases of density modulatlont where density modulation 
is different in lateral extent by layers where high density over row 
centers falls to zero between rows. The adoption of a profile or shape 
precludes the modeling of gradual and partial plant extension into the 
inter-row region. The point of view taken here is Chat the canopy is 
still represented by Infinitely extended layers even if no plant components 
occupy parts of the layer. That is« a profile per se is not considered to 
be the feature causing the row effect in field reflectance bi- rather it 
is Che density modulation which may or may not result in a profile. 

In the following extension of the uniform canopy model* the density 
modulation will be the same for all layers so that a particular profile 
would not be evident to the eye as Illustrated in Figure 3. The use of 
the same density modulation for all layers simplifies the calculations 
but should still lead to the essential features of the row effect on 
canopy reflectance. 

3.2 ROW EFFECT LOGIC 

Let the row modulation be defined by 

M(6) - n(6)/n(IF0V) (15) 

where 6 « across-row strip displacement* 

n(6) ■ mean strip density of components* 
n(IFOV) “ mean value for the IFOV. 

Since the modulation should be periodic and since the average of 
n(6) across rows must converge to n(lFOV)* then* 

M(6) - M(P + 6) i 0* (16) 

P 

f M(6) d6 - 1* (17) 

O 

where P * the row spacing or period. 
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Let the five parameters* a* b* c* c' and k* be the IFOV mean values. 
Then the five parameters for strips required for row structure must be 
simply the IFOV means multiplied by the modulation* M(5)* since all para- 
me;.ersvary in direct proportion to component density. 

Now* using the same differential equations as before but with the 
five parameters required for row structure* one obtains* 


dE(+d)/dz - -M(6) aE(4d) + M(6) bE(-d) 
+ M(6) cE(s)* 

(18) 

dE(-d)/dz - M(6) aE(-d) - M(6) E(+d) 
- M(6) c'E(s)* 

(19) 

dE(s)/dz ■ M(6) kE(s)* 

(20) 


for each strip at level z in the canopy displaced from the row center 
by distance* 6. 

The relations (18)* (19)* and (20) are ;o be solved and applied as in 
relation (11). Then the lateral average over all displacements* 6* must 
be calculated to find the average radiance in the direction of view. In 
relation (11) » the values of u* v* and w' are the IFOV means which must 
also be multlpled by the density modulation to account for density 
variations. Thus, relation (11) will become, 

P 

L - (1/P) y L(fi) dfi, (21) 

o 


where 

L(6) 
The field 


■ ^ Py^^v’ M(6)[uE(+d) + vE(-d) + w'E(8,6)l dz/n 


+ -h. 6)(E(-d* -h) + E(s, 5)1 p(soil)/ir 

reflectance is found as bcf'”*e. 
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The details of the solution to these equations and the determination 
of the average radiance are treated In Appendix A. In order to solve these 
equations! an Important assumption Is made that the upward and downward 
diffuse flux averaged along a strip is approximately uniform over all strip 
displacements. For small variations due to density modulation! this 
assumption is not different from the tacit assumption made in the uniform 
canopy model that the presence of random clumps will not result in signifi> 
cantly non-uniform lateral distributions of diffuse flux. The diffuse 
nature of this flux would tend to make the lateral distribution uniform in 
spite of local variations in concentration. 

However! some important agricultural canopies! such as! soybeans! have 
density modulations with large values along row centers falling to zero 
between rows. Unless some modification is made in the logic! the upward 
diffuse flux will not even be approximately laterally uniform. Visualize! 
for example, an extreme case where a strip of dense vegetation occupies 
the row center and several meters of bare soil occupy the inter -row region. 
The upward diffuse flux at soil level below the dense vegetation is very 
small but the upward reflected diffuse flux over the exposed and fully 
illuminated soil between rows is very large. The upward diffuse flux 
cannot be considered uniform by any stretch of the imagination. But the 
overwhelming majority of the flux reflected from inter-row bare soil does 
not participate in the canopy scattering process which is controlled by 
relations (18) and (19). Sunlight is Incident without scattering and 
is reflected largely towards outer space without coming close to any canopy 
components. Such flux is purely the result of a sun-soil interaction and 
would occur without any canopy at all. 

In order to remove the excess, purely sun-soil up^fard diffuse flux 
from the lower boundary condition, an estimation formula is used to determine 
the amount of this excess and the excess is removed from the upward diffuse 
flux at the lower boundary. The remaining diffuse flux which does parti- 
cipate in the canopy scattering is uniform relative to the total diffuse 



flux Involved in reflection from the IFOV. The contribution of the excess 
to IFOV radiance Is accounted for in the second term of relation (21) . 

With the above modificatloni the row effect calculation should apply to 
large, as well as small, variations In density. 
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4.0 

COMPARISON OF RESULTS WITH THE VERHOEF-BUNNIK MODEL 

The Verhoef-Bunnik rectangular prism profile model for calculating 
the row effect on crop reflectances is a special case of large density 
modulation where M(6) varies between some constant value over the vegeta- 
tion to zero over the bare inter-row soil. They calculated the variation 
of field reflectance with direction of view for a hypothetical field consisting 
only of green wheat leaves. Fortunately, they also published the necessary 
parameters by which a compariaor calculation can be made using this model. 

The results of a comparison calculation are shown in Figures 4, 5. 
and 6. The figures are polar plots of reflectance as a function of view 
azimuth. The row direction is north-south on the plot. In these figures, 
the rectangle points are the predictions of the uniform canopy model where 
vegetation is distributed randomly over the field. The "x" points are the 
predictions of the Verhoef-Bunnik rectangular prism nndel and the solid 
dots are the predictions of the new model. The closest agreement occurs 
when the view direction is close to the row direction. The variation in 
reflectance compares fairly well in other directions also except on the 
sunlit side of the canopy for 550 nm and 670 nm. Calculations for other 
sun directions yield similar results with some disagreement on the sunlit 
side. Because of the diversity in the methods of calculation in these 
t%ro models, the cause of the disagreement on the sunlit side could not be 
traced. The validity of predictions will depend upon field tests. 

Both models predict similar qualitative trends, however. In particular, 
both predict that row structure is most evident in the chlorophyll absorp- 
tion band. This variation is primarily due to the large contrast between 
low green leaf reflectance and high soil reflectance that very often 
occurs at 670 nm. 
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POLM SUN - 30* 

POLAR VIEW • «S* 

RECTAN6LE • UNIFORN 
X • VERHOEF AND lUNNIK ROM 
SOLID DOT • SUITS ROM 



FIGURE 4. (XmPARISON OF PREDICTIONS FOR 550 nn. 
The reflectance as a function of view azinuth 
shows the comparison between the Verhoef-Bunnik 
model and this row model for the rectangular 
prism profile ffheat leaf canopy. 
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FIGURE 5. COMPARISON OF PREDICTIONS FOR 670 nm. 
The reflectance as a function of view azimuth 
shows the comparison between the Verhoef-Bunnik 
model and this row model for the rectangular 
prism profile wheat leaf canopy. 
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FIG'" J 6. COMPARISOT OF PREDICTIONS FOR 870 nm. 
the reflectance ea a function of view arlmith 
8ho%fS the comparison between the Verhoef-Bunnlk 
model and this row model for the rectangular 
prism profile wheat leaf canopy. 
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Both predict very little effect due to row structure at 870 nin. The 
minor variation is due to the low soil to leaf contrast, as well as the 
high value of diffuse flux which tends to produce uniform illumination 
in spite of row structure. 
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5.0 

ROW MODEL PREDICTIONS FOR WHEAT 

One of the purposes of canopy reflectance modeling is to assess the 
effect'’ • canopy structure upon received signals* The crop which is 
considered here to assess row structure effects is wheat. The spectral 
properties of plant components and soil were those for Kansas wheat and 
soil obtained under previous NASA contracts. The model equivalent compo- 
nent sizes were obtained from Ray Jackson's field data for wheat as a 
function of Feekes scale. 

Two wheat stages were modeled: Feekes 5 and Feekes 8. The row 

modulation was taken to be the rectangular prism modulation of Verhoef 
and Bunnik which might be suitable at Feekes 5 but probrbly is not accurate 
for Feekes 8. Figures 7, 8, and 9 show polar plots of reflectance for 
three band centers — 550, 650, and 750 nm. Row direction is north-south 
in the plot and the direction of view is the nadir in all cases. Because 
of the symmetry due to the nadir view, only one sun azimuth quadrant for 
each band center is necessary to illustrate all of the important variations. 
Along with solar azimuthal variations shown on the polar plot, three 
different polar sun angles were used. The solid polar plot is for a 
25^ polar angle, the long dash plot is for a 45° polar angle, and the 
short dash is for a 60° polar angle. The scale for the 750 nm band is 
different from the scale for the 550 and 650 nm bands. 

Figure 7 shows the results for Feekes 5 wheat. One may observe for 
variable sun angle similar qualitative features as in the Verhoef- 
Bunnik calculation for variable view angles. The greatest effect is in 
the 650 nm band and the effect becomes more significant as the polar sun 
angle increases while the infrared 750 nm band is only moderately affected. 
One can see that the infrared to red ratio, which is often used as a 
crop vigor measure, will be significantly altered merely by sun-to-row 
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FIGURE 7. RECTANGULAR ROW WHEAT, FBEKES 5. 
Rows are north-south and the view direction 
Is the nadir. Reflectance In three bands Is 
shown In a polar plot as a function of sun 
azimuth for three polar sun angles: solar 

line “ 25”, long dash « 45®, short dash - 65®. 
No sky light. 
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FIGURE 8. RECTANGULAR ROW WHEAT, FEEKES 8. 
Rows are no^th-south and the view direction 
is the nadir. Reflectance in three bands is 
shown in a polar plot as a function of sun 
azimuth for three polar sun angles: solid 

line - 25®, long dash - 45°, short dash ■ 65°. 
No sky light. 
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FIGURE 9. MODIFIED ROW WHEAT. FEEKES 8. 

Rows are north*south and the view direction 
is the nadir. Reflectance in three bands is 
shown in a polar plot as a function of sun 
azimuth for three polar sun angles: solid 

line ■' 25®* long dash * 45®* short dash * 65®. 
No sky light. 
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angle conditions. These calculations are for direct sunlight alone. 

The addition of skylight will tend to reduce the extreme variations for 
the setting sun. 

The case for Feekes 8 wheat for sunlight alone is shown in Figure 8 
for an extreme rectangular prism type of row structure which would almost 
surely require artificial restraints to hold the vegetation away from the 
inter-row region. The row effect is also extreme. Figure 9 shows the 
same wheat where the row structure allows only 5X of the peak on-row 
concentration to appear at mid-row. Notice that the row effect is still 
significant but is much more subdued. It does not take much more vegetation 
in the inter-row region to reduce the row effect to negligible proportions. 

The impact of row direction on Landsat signals from the latter field 
was estimated for a 45° sun angle and a nominal amount of path radiance. 

The resulting MSS7/MSS5 ratio and Greenness measures are shown in Table I 
for sun down-row and sun across-row directions. 


TABLE I. 

ESTIMATED LANDSAT RESULTS 


ACROSS-ROW 

DOWN-ROW 

MSS7/MSS5 

2.0 

1.33 

GREENNESS 

47.5 

42.1 


The down-row direction gives an indication of a much less vigorous 
field. An underestimation of crop vigor and biomass could result purely 
from a chance row-sun relation. However, the cross-row direction does not 
lead to a serious overestimation. The reflectance for the cross-row 
direction is not greatly different from that of the uniform canopy. 
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APPENDIX A 

CALCULATION DETAILS FOR ROW CANOPY 

The first step in calculating the canopy reflectance is to solve 
relations (18) t (19)* and (20) to determine the diffuse and specular 
flux at any level, z, within the canopy. Relation (20) for specular 
flux is solved first and substituted into relations (18) and (19). 

SOLVING FOR SPECULAR FLUX 

In order to solve for specular flux, one must recognize that the 
density modulation for slant rays through the canopy will also be a 
function of z. Consider a coordinate system with z > 0 at the top of 
the canopy, x ■ 0 at a position, 6, in the direction across rows, and the 
y axis along the row direction as shown in Figure lA. A unit vector 
directed at the sun, s’, will be given by, 

s (sin 6 sin ip , sin 6 cos , cos 6 ) (lA) 

s s s s s 

where 0^ ■* polar sun angle, 

^ * sun azimuth relative to row direction or y axis 

(positive clockwise). 

The density modulation along the incoming ray of sunlight will be 
M(6 + x^) but from relation (lA) the lateral distance, x^, along a ray 
is given by, 

X ■ (z - z) tan 6 sin 

r r s s (2A) 


or 


X ■ a (z 
r 8 r 


z) 


where a “ tan 0 sin \b , 

S S ^8 

x^ * X coordinate along the ray, 
z^ « z coordinate along the ray. 
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Consequently, the density modulation along a ray is M(6 ^ ~ >))• 

Relation (20) may be written as, 

dP.(«. *)/d* ■ kM(« + a (* - t)) p_(fi, z) (3A) 

w 9 R 9 t 9 O 

where p^(5> z) * probability of line of sight to level z in 
the direction, T, 

■ specular irradiance at the top of the canopy. 

The Integration of (3A) is then, 

o 

Pg(6, z) - exp(ky^ M(6 + ***r) 

z 

and E(s) - E p (6, z). (5A) 

O 9 

Numerical methods will be needed for the integration indicated in 
relation (4A) since any density modulation may be considered. 

It should also be evident that the probability of line of sight, 
p^(6, z) , in the direction of view will be obtained in like manner but 
with the unit vector pointing towards the viewer with a ~ tan 0 sin 0 / 

V V V 

and with extinction coefficient, K. 

The noxt step is to substitute E(s) into relations (18) and (19) 
and solve. The assumption is made that E(4d) and E(-d), because of their 
diffuse nature, are approximately laterally uniform even with density 
modulation. The diffuse flux is controlled by the lateral average of 
diffuse extinction, diffuse scattering, and specular flux scattering. 

Thus relations (18) and (19) become, 

dE(-Kl)/dz - -aE(-W) + bE(-d) + cE M(6) p (6. z) . (6A) 

o s 

dE(-d)/dz - aE( H) - bE(4d)-c"E M(6) p (6, z) , (7A) 

o s 

where the bar indicates the lateral average over 6 . 
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Again, nuDwrlcal methods are required for Che indicated integration for 
Che average. 

The solution to (6A) and (7A) are now found by standard procedures. 

The diffuse flux variables are separated by use of differential operators 
and the resulting equations are solved by finding the solution to the 
homogeneous equations and adding a special solution using Lagrange's 
method of variation of parameters (also known as variation of conatants). 

In solving for the upward and downward welling diffuse flux, the 
lower boundary condition normally requires that, 

E(-hl, -h) - p(8oil)(E(-d, -h) + E p (6, ^h)). (8A) 

o s 

However, part of the contribution from E p (6, -h) represents excess 

o s 

sun- 80 ll interaction and does not contribute to the ‘upward diffuse flux 
within the canopy. An estimation of the amount of the excess is made by» 

E (sun) ■ p(soil) E Ip (-h, 6)t(6) - p (-hp uniform) T(uniform)J (9A) 

X 0^8 S ^ 

where p (-h* 6) • fraction of direct sunlight on the soil at 

lateral position^ 5 

t(5) ■ diffuse transmittance of the canopy if the entire 
canopy had the modulation density above positionp 
6* with all components black and opaque. 

p (-hg uniform) ■ fraction of sunlight on the soil for a uniform 
^ canopy with IFOV average. 

T (uniform) • the diffuse transmittance of the uniform canopy 
with all components black and opaque. 

The bar indicates the lateral average as before. 

One may observe that the excess is zero for a uniform canopy. The 
two terms in relation (9A) represent the estimated amount of direct 
sunlight which both reaches the soil and also is able to escape as diffuse 
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flux without participation in canopy acattaring. Tha aacond tarn aatiatatea 
the noroMil amount of aacapt and tha firat tarm aatimataa tha amount aacaping 
due to row structura. 


The adjuatad lower boundary condition becomes • 

E(-Hi, -h) - p(soil) [E(-d, -h) + E^ p^(«. -h)] - E^(sun) 


(lOA) 


The final step is the calculation of tha lateral average of observed 
radiance using relation (21). Reversing the order of integration! one 
obtains. 

L ■ j M(6) p (6, z)[uE(4d) + vE(-d)) dz/it 
-h ^ 


/ '' 


+ / w'E M(6) p (6, *) p (6, z) dz/ir 
~h ” 


+ pfsolDfp (6. -h) E(-d, -h) + E p (6, -h) p ( 6 , “h)|/ir (llA) 

L V 0 8 V J 


Numerical methods must be used since M(6) may have any periodic variation 
desired. 

For skylight alone, E ■■ 0, and the lower boundary condition is. 


E(-d, o) - E(sky). (12A) 

But some fractional part of the diffuse flux at soil level is only 
sky-soil interaction due to row structure. The remaining diffuse flux is 
spectrally modified by canopy scattering. The excess fraction of spectrally 
unmodified skylight at soil level is estimated by, 

E^(aky) ■ ^ 1^(6) - T^(uniform)j E(aky) (13A) 
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Thus I the rediencc contribution of soil illunineted by diffuee flux is* 

L(aoil) - p(«oil)|E(-d, -h)|l - [t^(6) - T2(uniforai)J J 

^ , (14A) 

+ - T2(uniform)jE(eky)j^Tr 

Both of the estimation formulas for direct and diffuse excess due 
to row structure are not mathematically derivable but are chosen on the 
basis of reasonableness in order to extend the scope of the model to 
Include extreme row structure circumstances where diffuse flux flow is 
actually no longer a one dimensional flow as assumed. 

Since numerical integration is required in a number of steps, 
computer programming of the calculation is necessary All numerical 
integrations that are independent of wavelength should be performed 
first and stored in computer memory as look-up tables to be used for 
the spectrally dependent calculations. 

The program is complex but it will easily fit on a personal computer 
with 12K of available RAM and run in BASIC. 


